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Abstract 

Local isotropy, or the statistical isotropy of small scales, is one of 
the basic assumptions underlying Kolmogorov's theory of universality of 
small-scale turbulent motion. The literature is replete with studies pur- 
porting to examine its validity and limitations. While, until the mid- 
seventies or so, local isotropy was accepted as a plausible approximation 
at high enough Reynolds numbers, various empirical observations that 
have accumulated since then suggest that local isotropy may not obtain 
at any Reynolds number. This throws doubt on the existence of universal 
aspects of turbulence. Part of the problem in refining this loose statement 
is the absence until now of serious efforts to separate the isotropic com- 
ponent of any statistical object from its anisotropic components. These 
notes examine in some detail the isotropic and anisotropic contributions 
to structure functions by considering their S0(3) decomposition. After an 
initial discussion of the status of local isotropy (section 1) and the theo- 
retical background for the SO (3) decomposition (section 2), we provide an 
account of the experimental data (section 3) and their analysis (sections 
4, 5 and 6). Viewed in terms of the relative importance of the isotropic 
part to the anisotropic parts of structure functions, the basic conclusion 
is that the isotropic part dominates the small scales at least up to order 
6. This follows from the fact that, at least up to that order, there exists 
a hierarchy of increasingly larger power-law exponents, corresponding to 
increasingly higher-order anisotropic sectors of the SO(3) decomposition. 
The numerical values of the exponents deduced from experiment suggest 
that the anisotropic parts in each order roll ofi' less sharply than previously 
thought by dimensional considerations, but they do so nevertheless. 

1 Introduction 

Local isotropy, or the isotropy of small scales of turbulent motion, is one of the 
assumptions at the core of the belief that small scales attain some semblance 
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Figure 1: The magnitude of the skewness of the temperature derivative in 
turbulent shear flows. The figure is taken from Ref. |2^. The Taylor mi- 
croscale Reynolds number R\ is proportional to the square root of the large 
scale Reynolds number. It is defined more precisely in section 2. 



of universality This is a statistical concept, and is not necessarily opposed 
to the idea of structured geometry of small scales |9| . The important question 
about local isotropy is not whether the small scales are strictly isotropic, but the 
degree to which the notion becomes a better approximation as the scales become 
smaller |2^ . Aside from the generic requirement that the flow Reynolds number 
be high (so that small scales as a distinct range may exist independent of the 
large scale), the question of how or if local isotropy becomes a good working 
approximation depends on the nature of large-scale anisotropy, on whether or 
not there are other body forces, on the nearness to physical boundaries, etc. All 
of this was well understood by the time of publication of Monin & Yaglom . 
The general consensus at the time seems to have been that small scales indeed 
attain isotropy far from the boundary at high enough Reynolds numbers, at 
least when one considered second-order quantities. 

The situation changed perceptibly when the small scales of scalar fluctua- 
tions were experimentally found to be anisotropic in many types of shear flows 
even at the highest Reynolds numbers of measurement. For a summary, see 
Ref. Since that time, various other pieces of evidence are slowly accumu- 
lating to suggest that small-scale velocity is also anisotropic; see, for example, 
Ref. The claim is that the previously held belief — which, to some de- 

gree, was comforting — loomed large only because we had not explored the right 
statistical quantities. The notion that anisotropy persists at all scales at all 
Reynolds numbers (though manifested only in certain statistical parameters) 
puts a strong damper on any theory that purports to consider small-scale tur- 
bulence as a universal object. This is somewhat of an impasse. 

Since the problems began with passive scalars, we might as well consider the 
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evidence in that case a little more closely. The evidence, collected over many 
years by many people, is reproduced in Fig. |l| from Ref. Q. The argument is 
that if the temperature fluctuation 9 is locally isotropic, its derivative dO/dx, 
being a small-scale quantity by construction, should be isotropic. Taking reflec- 
tion symmetry as part of isotropy, we should have {{dO/dx)^) = —{{d9/dx)^) 
for all odd values of n. This means that all the odd moments must be zero. In 
particular, {{dd/dx)^) = 0, or 'small' in practice. But small compared to what? 
The standard thing to do is to normalize {{dO/dx)'^) by {{89 / dx)'^)^^'^ , or to 
examine the behavior of the skewness S of dO/dx. These are the data shown 
in Fig. |l|. The data suggest, despite some large scatter, that the skewness re- 
mains to be of the order unity even at the largest Reynolds numbers for which 
measurements are available. 

There are two points to be made. The flrst is that the Taylor microscale 
Reynolds number R\ may not be the right parameter against which to plot S. 
As articulated by Hill Q , the reasons are that the velocity that appears in R\ 
is a large-scale quantity and the length scale A cannot be defined independent 
of the large scale velocity. (In practice, the definition of A may not require the 
large scale, see , though this is not a result that can be shown formally; but 
the comment on the velocity remains valid in any case.) We believe, however, 
that the situation will probably not change qualitatively even if we adopted a 
different abscissae for Fig. |l|. The second point is that, if one used {{39 / dx)'^)'^/'^ 
instead of {{89 /dx)'^)^^'^ for normalizing {{89/8xY), the normalized quantity 
will vanish with Reynolds number, roughly according to some power of the 
Reynolds number. Why is it not legitimate to compare the third moment to the 
next high-order even moment, instead of the neighboring low-order moment, or, 
perhaps to the geometric mean of the lower and higher order even moments? In 
either case, we will have a quantity that diminishes with the Reynolds number. 

These are not elegant arguments, and seem like desperate efforts made to 
save the situation at any cost. There is a further argument to be made, however. 
Leaving aside the technicalities for a moment, let us suppose that, for a given 
statistical quantity to be measured, there are at all scales an isotropic part and 
an anisotropic part. We might now ask whether the ratio of the anisotropic 
part to the isotropic part vanishes as the scale size becomes smaller. We are no 
longer asking if the third moment vanishes with respect to (a suitable power of) 
some even order moment, but if, within a moment of a given order, the isotropic 
part eventually dominates the anisotropic part as the scale size vanishes. If this 
is indeed so, it allows us to say that, while non-universal anisotropic parts may 
always be present, the isotropic part dominates at small enough scales. If, of 
course, the isotropic part contributes exactly zero to the statistical quantity be- 
ing considered, the anisotropic part will no doubt prevail at any flnite Reynolds 
number, but this is not in contradiction to the statement just made. This is 
a richer point of view to take, potentially less inhibiting for the development 
of a sensible theory of turbulence. In such a picture, the universal theory that 
may emerge holds for the isotropic part alone, but, for it to be applicable to 
other types of turbulence, the anisotropic 'correction' has to be added in some 
suitable way. 
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The questions of interest, then, are obvious: What is a good way to de- 
compose any statistical object of choice into isotropic and anisotropic parts? 
How do these two parts vary with scale size, relative to each other? A plausi- 
ble method for answering these questions was proposed in Ref. Q by using the 
S0(3) decomposition of tensorial objects usually considered in turbulence. One 
can argue as to whether this is the best perspective, but there is no doubt that it 
provides one framework within which our questions can be posed and answered. 
An experimental assessment of this issue is the broad topic of these notes. This 
is preceded by a detailed description of the theoretical issues involved. 

The notes were part of the lectures presented by KRS at Les Houches, and 
form part of the Ph.D. work of SK. We note in advance that the literature cited 
is limited to a few key articles. We thank Itamar Procaccia and Victor L'vov 
for introducing us to the subject discussed here, and Christopher White and 
Brindesh Dhruva for their help in acquiring some of the data. KRS thanks the 
organizers of the School for the invitation to deliver the lectures. 

2 Theoretical tools 

2.1 The method of SO(3) decomposition 

A familiar example of decomposition into the irreducible representations of the 
S0(3) symmetry group is the solution of the Laplace equation in spherical co- 
ordinates 



where "0 is a scalar function defined over the sphere of radius r. This equation 
is linear, homogeneous and isotropic. The solutions are separable, and may be 
written as products of functions of r, and (p. The solutions form a linear space, 
a possible basis for which are derived from the spherical harmonics YimiO, ip) 



where the index I = 0, 1,2, ...,cxo denotes the degree of the harmonic polyno- 
mial. The angular dependence {9, (p) is equivalent to the unit vector f , or the 
unit sphere. In the angular dependence, the index m denotes the elements of 
the orthonormal basis space that span each harmonic polynomial of degree I. 
For each I there are 21 + 1 elements indexed by m = — j, —j + 1, . . . , j which 
are the Yi„i. Equation (Q) is a useful basis space of solutions to the Laplace 
equation because the -tpim (f ) are each solutions of the Laplace equation and are 
orthonormal for different l,m. Each basis element has a definite behavior under 
rotations, that is, the action on it by an element of S0(3) preserves the indices. 
In other words, the rotated basis element will also be indexed by the same 1,171. 
The general solution of (0) is given by 



VV(r) = 0, 



(1) 



(2) 




(3) 



l,7n 
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The coefficients aim are obtained from the boundary conditions on ^/;(r) and, 
for a particular solution, any of them may be zero. In the theory of group 
representations, V'im(r) is said to be the 21 + 1 dimensional irreducible represen- 
tation of the group of all rotations, the S0(3) symmetry group, in the space of 
scalar functions over the sphere. The space of scalar functions is called a 'car- 
rier' or 'target' space for the S0(3) group representations. For each irreducible 
representation indexed by I = O...00 there are 21 + 1 components indexed by 
m = —I... + L A representation of a symmetry group is said to be irreducible if 
it does not contain any subspaces that are invariant under the transformation 
associated with that symmetry. 

Here we make a distinction between the behaviors under rotation of the equa- 
tion (|l|) and of its solutions ip{r). The statement that the equation is isotropic 
means that it will hold true in a rotated frame with the Laplacian operator V and 
the function ^(f) properly defined in the new coordinate system, V ip{^) = 0- 
However, the function V'(r) may contain both isotropic and anisotropic com- 
ponents. All symmetry groups possess a one-dimensional representation in the 
carrier space, which is invariant under the transformations of that group. For 
the representations of the S0(3) symmetry group in the space of scalar functions 
over the unit sphere, the one-dimensional representation, indexed by Z = 0, is 
Koo = constant] it will clearly remain unchanged under proper rotations. This 
is the isotropic representation. The higher-dimensional representations Yim (di- 
mension 3, 5, 7, . . . corresponding to / = 1, 2, 3, . . .) have functional forms in 9 
and (p which are altered under rotation even though the degree (1), and hence 
the dimensionality, of the representation is preserved. These are the anisotropic 
representations. We examine this point in more detail with further discussion 
of the simple example of the scalar functions ip. 

For a particular rotation A in Euclidean space which tells us how to rotate 
a vector x into a new coordinate system where it is denoted by x, 

x"=A'lx'^, (4) 

we can define an operator 0\ which tells us how to rotate the function ipir) to 
V'(f). If the function V'(r) is written in terms of a linear combination of its basis 
elements as in Eq. (|^), then the rotation operation is written as 

V^(f) - OA^(r) 

= OA^a/mV/m(r) 

= ^OAa/„'i/'/m(r)- (5) 

The Oa for each I representation is a {21 -I- 1) x {21 + 1) matrix denoted by 
^m' m(A)- transformation is written as 

I 

OAi^Uv) = Di'Xmim^'ir)- (6) 

m' — — l 
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Thus, when a rotation of the function ip(r) into a new coordinate system is to 
be performed, the rotation matrices dI^^,{A) are aU that is needed in order 
to transform each of the irreducible representations of the SO (3) group that 
form the basis functions of i^{r). As can be seen from Eq. (^), the j — 
(one-dimensional) irreducible representation is the one which is invariant to all 
rotations in that its functional form is always preserved. The _D^*'-'(A) matrix 
is merely a number, independent of A, multiplying the j — representation. 
The j = irreducible representation is the isotropic component of the SO (3) 
decomposition. For all higher-order j's the rotation preserves the dimension 
of the transformed component, i.e. the indices j and m are retained, but the 
functional form is altered. The D(^^(A) for j > are true (2j -I- 1) x (2j -I- 1) 
matrices that mix the various m contributions of the original, unrotated basis. 

The above simple example involved the case of a scalar function ip that 
depended on a single unit vector f . The method of S0(3) decomposition may 
be carried over to more complicated objects. The theory is given in detail in 
. In general we can imagine an n*'* order tensor function which depends on p 
unit vectors j'aias - "" (f i, f2, . . . fp) and is the solution to an isotropic equation. 
Then, the rules of SO (3) decomposition carry over in the following manner. The 
rotation operator 0\ is now defined through the relation 

_^ ^ _^ ^ Q^j,aia....a„ (f fa, . . . , fp) 

= A''pIA''^...A$:T''^^^---^'-{A-'H,A-'I^,...,A-%). (7) 

The tensor function T may be written in terms of the irreducible representations 
of the S0(3) symmetry group. The basis elements are denoted by i?"^"^ '"" (f ) 
where the j index plays the same role as the I index in the example of the 
scalar function. The basis elements are more complicated than the ipim because 
the space of functions T is the direct product of n Euclidean three-dimensional 
vector spaces (manifest in the indices ai...a„) with p infinite dimensional 
spaces of continuous single- variable functions over the unit sphere (the f i . . . fp). 
If the constituent spaces are also written in the SO (3) decomposition, the rules 
of angular momentum addition, familiar from quantum mechanics, may be used 
in taking the direct product. The three-dimensional Euclidean vector space is a 
J = 1 space while each of the infinite-dimensional spaces functions over the unit 
sphere is the sum of the j = 0, j = 1, . . . ,j — > oo irreducible representations of 
S0(3) (recall the above example of V'(f )) with each j representation appearing 
once. 

In tensor product notation, the product space of two vector spaces Vi and 
V2 is denoted by Vi (E) V2- In the case of vector spaces in the SO (3) notation 
the spaces are named uniquely by their index j. Using tensor product notation, 
the n three-dimensional Euclidean spaces form a space 1 (g) 1 (g) ... 1 (n times) . 
The tensor notation indicating a linear sum of tensors with SO (3) representation 
indices ji and j2 is ji ® j2- Using this notation, each of the infinite dimensional 
spaces is ® 1 © 2 © . . .. The direct product of p of these is (0 1 2 ® . . .) (g) 
(0 © 1 © 2 ...)©...© (0 ® 1 © 2 e .. .) (p times). Thus, the final product space 
for Bjm (f ) is written in tensor-product notation as 1 (g 1 (g ... © 1 {n times) 
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(g)(0©l©2©. . .)(8>(0©1©2©. . .) ■ . • (0©1©2©. . .) (p times). Now, the tensor 
product of spaces ji and j2 will contain new spaces whose S0(3) indices are given 
in the following manner. We recall the rules of angular momentum addition 
familiar from quantum mechanics. The total angular momentum of an SO (3) 
representation space is given by its index j. The rules of angular momentum 
addition dictate how the product of two spaces of angular momentum (SO (3) 
index) ji and j2 may be added and the possible j indices of the resulting spaces, 

ji <S> j2 = |ji - j2 1 © . . . © (ji + j2). (8) 

Equation (||) says that the direct product of two spaces each belonging to a 
particular j will generate a sum of new spaces with only those j indices allowed 
by the rule. The operations (Ei and © are distributive like the corresponding 
arithmetic operators. For example 

10 1 = 0© 1 ©2 
l©l(g)l = (0©1©2)(8)1 

= (0®1)©(1®1)©(2©1). (9) 

If we now apply the angular momentum addition rules (^) to each of the product 
terms on the right hand side of Eq. , we get 

l©l(g)l = 1©0©1©2©1©2©3 

= 0© (1 X 3) © (2 X 2) ©3. (10) 

When taking the product of more than two spaces, there will be several ways 
to arrive at a particular j in the final product space. We see this in the final 
rearrangement of the right hand side of Eq. (p^, where there is only one repre- 
sentation each of j = and j — 3, but three of j = 1 and two of j — 2. Going 
back to our more complicated target space for the basis elements B"^^"^ ' "" (f ), 
which is l©l(g). . .©1 (n times) ®(0©1©2©. . .)©(0©1©2©. ..)... (0©1©2©. . .) 
(p times), we expect that there may be many ways to obtain a particular j in 
the final product space. To indicate this, for each basis element we associate 
a further index q giving Bqjm{^)- These are important when we start to actu- 
ally calculate the basis elements. The machinery used is the Clebsch-Gordan 
method well-known in quantum mechanics and the reader is referred to |^ for 
the details. The tensor T may also depend on vectors ri with magnitudes but 
for each rotation operation Oa these may be treated as parameters included in 
the weight associated with each basis element. The T are represented as 

r"^""-""(ri,r2,...,rp) = ^ a,,„.(ri, ra, . . . , r^) 

XB,",C -""(fl,f2,...,fp). (11) 

In the next two sections, we move away from the very general formalism of 
S0(3) decomposition discussed above, and apply it to the specific case of the 
statistical tensor quantities in fluid turbulence and their dynamical equations 
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obtained from the Navier-Stokes equations. We are motivated by the observa- 
tion that the S0(3) decomposition when apphed in the correct manner would 
allow (a) separation of the scaling variable r from the angular dependence, and 
(6) the separation of the isotropic from the anisotropic parts of the tensor. 

2.2 Foliation of the structure function into j-sectors 

We would like to use the formalism of the SO (3) representation, as presented in 
general terms in the previous section, to study the structure function tensor in 
turbulence theory. The nth order structure function is defined by an ensemble 
average of the nth moment of the difference of velocity components across scales 
Tp. It is in fact a tensor function over p unit spheres. We define the velocity 
difference as 

«;"'(rk) = u"Kx + rk)-w"'(x), (12) 

where the ai denotes the component of the velocity vector in the direction j in a 
defined coordinate system, and the subscript k on rt denotes a particular choice 
of the vector. The nth order structure function for p such choices of vector scale 
r is then 

^aia....a„ (j.^^ j.^^ . . . , rp) = {t,)w^- (rj) . . . x i/;"" (rO), (13) 

where the subscripts ■ ■ denote any of the vectors from ri to Fp. 

In this section we review briefly the reason why the structure functions 
may be written as a linear combination of the irreducible representations of 
the S0(3) symmetry group. For details the reader is again referred to Q. 
The dynamical equation for the structure function may be derived from the 
Navier-Stokes equations. Similar to the solutions of the Laplace equation, the 
solution to the nth order structure function equation forms a linear space with 
the basis chosen to be the irreducible representations Bqjmi^) of the SO (3) group 
as shown in previous sections (see Eq. (pi]). Since the basis is orthonormal for 
different j, m, and the equation is isotropic, we obtain a heirarchy of dynamical 
equations each with terms of a given j and m. This demonstrates that the 
dynamical equations themselves do not mix the various j and m contributions, 
and that the dynamical equation for the isotropic part of the structure function 
is different from the dynamical equations for any of the higher j contributions. 
Further, it may be shown from solvability conditions on this set of equations 
that the scaling part of the function changes between j-sectors while remaining 
the same within a j-sector, for different m. This motivates the postulate that 
the different j sectors scale with exponents different from the isotropic ^n- We 
denote these by where subscript n indicates the rank of the tensor and 
superscript (j) indicates the number of the irreducible representation. We will 
leave out the superscript in the case of j — since then we recover the known 
isotropic scaling exponent C„. 
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2.3 The velocity structure functions 

In the previous section we provided a heuristic justification for the use of the 
S0(3) decomposition for this tensor using the fact that it is the solution of an 
isotropic equation. Now, we use the rules of SO (3) (see Eq. (^) decomposition 
to write it as 

5'"'"'-""(ri,r2, . . . ,rp) = ^ agj>„(ri, ra, . . . , r^) 

q.j,va 

X i?,",^::-""(fi,f2,...,fp). (14) 

In this form, as was demonstrated for the scalar function ^/i, the j = 
representation is the isotropic one. Its prefactor aq^^{rl^r2^ ■ ■ ■ ,rp) contains the 
scale dependence. For p > 1 there will be infinitely many ways to obtain a 
particular j for the product space as may be seen from the angular momentum 
addition rules in Eq. (^) and the direct-product representation. Therefore, q 
ranges from 1 to oo. This is a difficult hurdle computationally but, fortunately, 
most experimental and theoretical work deals with the dependence on a single 
vector, and p — I. In that case, and for rank n = 2, we obtain the second order 
structure function 5"'^(r) quite easily. 

2.3.1 The second-order structure function 

We consider the structure function of Eq. ( p^ ) for n = 2 and p = I. The tensor 
product space is 

l®l(g)(0©l©2®...) = (0©1©2)®(0©1©2©...) 

= (0«)0)©(0®1)©(0®2)©... 

© (10 0) ©(10 1) ©(1(8)2)©... 

© (2 0) ©(2 1) ©(2® 2)©.... (15) 

As demonstrated in ( [To| ) using the addition rule (^), there could be more than 
one way of obtaining a particular j representation is obtained in the product 
space. To count these for a given j, we had to add there the index q to the basis 



tensors (see for example Eq. (11)). We find for the second-order tensor over a 
single sphere that 

• j — has a total of 3 representations 

• j = I has a total of 7 representations 

• j > 1 has a total of 9 representations. 

The Clebsch-Gordon machinery tells us in addition about the symmetry (in 
the indices) and parity (in r) of each q contribution. Using this information, 
the terms of the basis may be constructed. As in the case of the scalar function, 
the -B^jCi orthonormal for different j,m and q. In practice, the Clebsch- 
Gordon method of constructing these objects is rather tedious and so we follow 
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the alternative method offered in [||. We make use of the Clebsch-Gordon 
methods only to obtain the number of representations in each j, and its parity 
and symmetry properties. Armed with this information, the B^fm ^ns-y be 
constructed by a more convenient means. The idea is to use the fact that we 
already know an orthonormal basis in the SO (3) representation for the scalar 
function over the sphere (Eq. (||)). We can now 'add indices' to this basis in 
a way to be described, while retaining the properties under rotation, in other 
words its j and m values. To add indices in a way that does not change the j,m 
values is to perform contraction with the objects 6°'^ ,r" ,€"^'' and the partial 
derivative operator da- This method automatically takes care of the j and m 
properties; and all that needs to be done now is to apply the above operators to 
obtain the different q terms with the right symmetry and parity properties. This 
gives us a complete set of basis elements which are different than the Clebsch- 
Gordan method. The orthonormality among different q for the same j and m is 
lost, but the orthonormality among different j and m is maintained because we 
start out with a basis Eq. (||) which already possesses these properties. However, 
different q elements are linearly independent and span a given j, m-sector. The 
details on how this is done using the rules from ||] is presented in appendix 
A. In what follows we simply write down the components calculated using that 
method. 

The second order structure function tensor is 

5"/3(r) = ((m"(x + r) - m"(x))(u'^(x + r) - m^(x))) (16) 
which we decompose using the SO (3) irreducible representations _B^j^(f) as 

5"'3(r) = 5j!fo(r) + 5;fi(r) + 5j!f,(r) + ... 

= a,,^{r)B;f^{a(3){r)- (17) 

A further constraint is provided by the incompressibility condition 

a„S'"^(r) 0. (18) 

Each j, m-sector must separately satisfy the incompressibility condition since 
taking the partial derivative preserves the rotation properties and does not mix 
the j, m sectors. Therefore, the incompressibility condition provides a constraint 
among the different q contributions within a given j, m-sector. For example, 
the j = contribution has three different representations {q = {1, 2, 3}), two of 
which are symmetric and one antisymmetric in the indices a, (3. By definition, 
the structure function is symmetric in the indices, therefore the antisymmetric 
contribution will not appear giving only two q contributions. The incompress- 
ibility constraint gives 

(19) 
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and 

5a(aioor'^^'5"'3 + 0200^'^^ ^) = 0. (20) 



Motivated by the expectation that the structure functions scale as powers of r 
in the inertial range, we assumed in Eq. ( pO| ) that the scale-dependent prefactor 
is of the form aq^^r^" where aq^^ is a flow-dependent constant. Equation ( pO| ) 
results in a relationship between aioo and a2oo giving the final form of the 
isotropic j — sector with just one unknown coefficient cq to be determined by 
the flow boundary conditions. That is, 



SfU^) = cor^^ [(2 + C2)<^"^ - C2—^\ . (21) 

Here, (2 ~ 0.69 is the known empirically known anomalous second-order scal- 
ing exponent. We would like to assume a similar scaling form for the prefactor 
O'qjmi'f) for j > 1. In such a formulation, there is a hierarchy of scaling expo- 
nents which we denote by ^ ^ C2 corresponding to the higher-order j sectors. 
Successive j's indicate increasing degrees of anisotropy. The following section 
provides a justification for such a classification of the scaling of the various sec- 
tors. It is the larger goal of this article to examine these aspects of the theory 
with the help of high- Reynolds-number data. 



2.4 Dimensional estimates for the lowest order anisotropic 
scaling exponents 

In this section we present dimensional considerations to determine the "classi- 
cal" values expected for (2^^ and (2^'' the spirit of Kolmogorov's 1941 theory 
(henceforth called K41). We work on the level of K41 to produce the value 
(^2°'* = 2/3. Ignoring intermittency corrections is justified to the lowest order 

because the differences between any two values (2"^ and (2 ^ for J j' &re 

considerably larger than the intermittency corrections to either of them. 

(2) 

It is easiest to produce a dimensional estimate for Q • One simply asserts 
[HI that the j = 2 contribution is the first one appearing in S'"^(r) due to the 
existence of a shear. Since the shear is a second rank tensor, it can appear 
linearly in the j = 2 contribution to S"^{r). We thus have for any m, —j < 
m<j, 

S';U^)^T"P^'^f{r,l). (22) 

Here T°'f^'>^ is a constant dimensionless tensor made of 6°"^, and bilinear 
contributions made of the three unit vectors p, rh, n, which form our coor- 
dinate system as defined in appendix A. The form of Eq. (|2^) means that the 
dimensional function f{r,t) stands for the response of the second-order struc- 
ture function to a small external shear. Within the K41 dimensional reasoning, 
this function in the inertial interval can be made only of the mean energy flux 
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per unit time and mass, e and r itself. The only combination of e and r that 
yields the right dimensions of the function / is e^/'^r^/^ . Therefore 

S'^l^{v) ^ T^0iS^e^/\^'^ . (23) 

We thus find a "classical K41" value of Q to be 4/3. 

To obtain the value of C,'^'' we cannot proceed in the same way. We need a 
contribution that is linear (rather than bilinear) in the unit vectors p, rh, n. 
We cannot construct a contribution that is linear in the shear and yet does not 
vanish due to the incompressibility constraint. Thus there is a fundamental 
difference between the j = 2 term and the j = 1 term. While the former 
can be understood as an inhomogeneous term linear in the forced shear, the 
j = 1 term, being more subtle, may perhaps be connected to a solution of some 
homogeneous equation well within the inertial interval. 

We therefore need to consider some quantity other than the shear which 
could contribute to the anisotropy. One invariance in the inviscid limit is as 
given by Kelvin's circulation theorem. In the so-called Clebsch representation, 
one writes the Euler equation in terms of one complex field a(r, i), see for 
example]!^. In the k-representation, the Fourier component of the velocity 
field u(k, t) is determined from a bilinear combination of the complex field 

u(k,<) = ^y"d3fcid3fc2*(ki,k2)a*(ki,i)a(k2,t) (24) 



*(ki,k2) = i(^ki+k2-(ki-k2)^j^^-^). (25) 



It was argued in that this representation reveals a local conserved integral 
of motion given by 

^ [ d^kka*{k,t)a{k,t). (26) 

Note that this conserved quantity is a vector, and so it cannot have a finite 
mean in an isotropic system. Consider a correction to the second order structure 
function due to a flux of the integral of motion tt. The dimensionality of tt is 
[tt] — [e^/'^/r^/'^], and therefore now the dimensionless factor is tt r^/'^/e^/'^. As 
we did with the shear in the case of the j = 2, we assume here the expandability 
of SttS at small values of the flux tt and find 

5;ii(r) -T^'^-^Tf-'r (27) 

where T"^"' is a constant dimensionless tensor that is linear in the unit vectors 
p, rh, n. We thus find the "classical K4r' value Cj^-* is 1. The derivation of the 
exponents ('2^'' and (^2^^ is given in p2[ . It is concluded that dimensional analysis 
predicts values of 2/3, 1 and 4/3 for C2"'-' with j = 0, 1 and 2 respectively. It is 
not known at present how to continue this line of argument for j > 2. 

For the higher-order structure functions 5'"i"2 - a„ j^j.^ where n > 2, similar 
dimensional analysis may be performed in order to find the K41 contribution to 
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scaling due to shear (corresponding to the j — 2 component). For each n, the 
lowest order correction to scaling that is linear in the shear is n/3 + 2/3. 

2.5 Summary 

The technique of SO (3) decomposition may be used in order to write the struc- 
ture function tensor in terms of its isotropic part, indexed by j = 0, and higher- 
order anisotropic parts, indexed by j > 0. The dynamical equation for the 
structure function of order n foliates into a set of equations each of different 
j, TO. This motivates the postulate that the different sectors scale differently. 
The theory also indicates that any scaling behavior would depend only the j 
index and be independent of the m within that sector. There is allowance for 
dependence on boundary conditions, specific kinds of forcing and so on, in the 
unknown coefhcients aqjm in the SO (3) expansion (Eq. (pi])). We first consider 
the second-order structure function in a detailed manner in the light of this 
group representation. The analysis may be implemented for structure func- 
tions of any order but the task becomes computationally and conceptually more 
difficult for n > 2. In section 6, we present a means of circumventing these 
problems. 

We have reviewed the theoretical estimates for the low-order anisotropic 
scaling exponents for the second-order structure function. The method may be 
carried over to the higher-order objects as well, and we have derived the lowest 
order shear-dependent scaling contributions to the n*'* order structure function 
to be 7-"/3+2/3_ Pqj. ^j^g second-order structure function, we have shown that the 
j = 2 contribution corresponds to a low-order shear-dependency. Again for the 
second-order structure function, we have presented a conserved quantity in the 
Clebsch-representation which provides the correct dimensional properties for the 
j = 1 scaling contribution to be r. The arguments used are purely dimensional 
in the K41 sense, and do not take account of anomalous scaling. We expect 
that real turbulent flows will exhibit anomalous scaling in the isotropic sector, 
and in every anisotropic sector in the SO (3) hierarchy. The issue of anomalous 
exponents in turbulence has now multiplied several-fold, to all the j sectors, in 
light of the apparent universality that this work suggests. Thus, the theoretically 
predicted anisotropic exponents are to be treated merely as estimates, especially 
for higher-order structure functions (n > 2), since the scaling anomalies are 
expected to increase with the order n, similar to the isotropic case. 

In subsequent sections we demonstrate the use of experimental data in order 
to test the predictions of the theory. In particular, we provide explicit calcu- 
lations of measurable tensor quantities and extract scaling corrections due to 
anisotropy. 

3 Some experimental considerations 



13 



3.1 Background 

Experimental studies of turbulent flows at very high Reynolds numbers are 
usually limited in the sense that one measures the velocity field at a single 
spatial point, or a few spatial points, as a function of time [^, and uses Taylor's 
hypothesis to identify velocity increments at different times with those across 
spatial length scales, r. The standard outputs of such single-point measurements 
are the longitudinal two-point differences of the Eulerian velocity field and their 
moments. In homogeneous and isotropic turbulence, these structure functions 
are observed to vary as power-laws in r, with scaling exponents Cn 0- 

Recent progress in measurements and in simulations has begun to offer in- 
formation about the tensorial nature of structure functions. Ideally, one would 
like to measure the tensorial nth order structure functions defined in Eqs. ( p^ 
and (p^. Such information should be useful in studying the anisotropic effects 
induced by all practical means of forcing. 

3.2 Relevance of the anisotropic contributions 

In analyzing experimental data the model of "homogeneous and isotropic small- 
scale" is universally adopted, but it is important to examine the relevance of 
this model for realistic flows. As we will demonstrate in the next chapter, 
our data we use exhibit anisotropy down to fairly small scales We have 

shown mathematically that keeping the tensorial information helps significantly 
in disentangling different scaling contributions to structure functions. In the 
light of the SO (3) representation of section 2, where it is shown that anisotropy 
might lead to different scaling exponents for different tensorial components, a 
careful study of the various contributions is needed. This is our goal in the rest 
of this article. 

3.3 The measurements 

In order to extract a particular j contribution and the associated scaling expo- 
nent, one would ideally like to possess the statistics of the velocity at all points 
in three-dimensional space. One could then extract the j contribution of partic- 
ular interest by multiplying the full structure function by the appropriate Bqjm 
and integrating over a sphere of radius r. Orthogonality of the basis functions 
ensures that only the j contribution survives the integration. One could then 
perform this procedure for various r and extract the scaling behaviors. 

The method just described was adopted successfully in using data from 
direct numerical simulations of channel flows. The Taylor microscale Reynolds 
number Rx for the simulations was about 70. This is not large enough for a clear 
inertial scaling range to exist. The authors of Ref. resort to extended self- 
similarity (ESS) in both the isotropic and the anisotropic sectors. Nonetheless, 
the results indicated a scaling exponent of about 4/3 in the j — 2 sector. While 
the experimental data are limited to a few points in space, and the integration 
over the sphere is not possible, we are able to attain very high Reynolds numbers 



14 



especially in the atmosphere under steady conditions {R\ « 10,000 — 20,000). 
Despite the advantage of extended scaling range, we are however faced with a 
true superposition of contributions from various j sectors with no simple way 
of disentangling them as was done with the numerical data. However, as we 
will show in section 6, we can make a judicious choice of the tensor components 
studied, and obtain access to the anisotropic contributions. 

The tensor structure of the velocity structure functions is lost in the compu- 
tation of the usual single-point single-component measurement of longitudinal 
and transverse objects. This is because the part of the expansion that is depen- 
dent on the angle 9 is hidden: the longitudinal and transverse components set 
the value of to a constant. The boundary-dependent prefactors now collapse 
to a single number as 

Sf£,ir) ^ ari^'\ (28) 

On the one hand, this is a simple expression. On the other, the angular de- 
pendence is completely lost and there is no formal difference from the isotropic 
object 

^;=fo«=cor«^ (29) 

In order to see the true tensor character of the structure function we need an 
angular variation of the scale separation r. In the atmospheric boundary layer 
which offers the highest Reynolds numbers available, the simplest configuration 
that would allow us to do this is a two hot-wire combination separated by 
distance A in the spanwise direction (y), orthogonal to the mean-wind. By 
Taylor's frozen flow hypothesis, such a set-up will provide two simultaneous one- 
dimensional cuts through the flow. Therefore, one can measure the correlation 
between the two probes, across a scale r that makes an angle 6 with the mean- 
wind direction. As r varies, the angle with respect to the mean-wind will also 
vary, giving some functional dependence on 9 for the coefficients in the SO (3) 
decomposition. A schematic of the experimental configuration is presented in 
Fig. |. 

A final consideration in measurements aimed at measuring anisotropic con- 
tributions is the homogeneity of the flow. The incompressibility condition may 
only be applied as a constraint on the structure function coefficients if the flow is 
homogeneous in the r direction. If we consider, for instance, probes separated 
in the shear (vertical) direction, r will have a component in the inhomoge- 
neous direction and incompressibility may not be used to constrain the terms 
in the S0(3) expansion, leaving too many free parameters. However, in some 
instances, it is important to consider quantities that are not constrained by the 
incompressibility condition. For example, the j = 1 contribution to the second 
order structure function will not appear in the experimental configuration de- 
scribed above. This is so because only one symmetric, even-parity, term exists 
in j — 1 by the Clebsch- Gordon rules, and this must vanish by the incompress- 
ibility constraint. Probe-separation in the shear direction would be needed to 
produce a non-zero contribution of the j = 1 term. In short, one can understand 
the nature of the best experimental configuration from carefully studying the 
tensor decomposition of the structure function. 
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_ The chosen coordinate system 



Figure 2: Schematic of the experimental set-up. Shown is the positioning of the 
single wire probes 1 and 2 with respect to the mean wind, and a visual expla- 
nation of how Taylor's hypothesis is employed. The coordinate system chosen 
is different from what is conventional in turbulence — where the mean-wind is 
taken to lie along the a;-axis. The pictured choice of coordinates simplifies the 
calculations involving the spherical harmonics. 
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mm 
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6 


4.1 


1.08 


1.1 


0.75 


15 


10,500 


10,000 


4 X 10^ 


35 


8.3 


2.30 


7.8 


0.45 


13 


19,500 


5.000 


4 X 10' 


0.11 


2.7 


0.47 


6.6 


0.47 


2.8 


900 


5,000 


8 X 10*^ 


0.27 


3.1 


0.48 


2.8 


0.6 


4.4 


1400 


5,000 


8 X 10^ 


0.54 


3.5 


0.5 


1.5 


0.7 


6.2 


2100 


5,000 


8 X 10^ 



Table 1: Data sets I (first line), II (second line) and III (third- fifth lines). The 

various symbols have the following meanings: U = local mean velocity, u' = 
root-mean-square velocity, (e) = energy dissipation rate obtained by the assump- 
tion of local isotropy and Taylor's hypothesis, rj and A are the Kolmogorov and 
Taylor length scales, respectively, the microscale Reynolds number Rx = u'X/u, 
and fs is the sampling frequency. 
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Figure 3: A typical compensated longitudinal third-order structure function. 
The inset indicates the region of isotropic inertial range scaling. (From Ref. fl.) 



We analyze measurements in atmospheric turbulence at various heights above 
the ground (data sets I, II and III) . Sets I and III were acquired from flow over 
a long fetch in the salt flats in Utah. The site of measurements was chosen to 
provide steady wind conditions. The surface of the desert was smooth and even: 
the water that floods the land during spring recedes uniformly and leaves the 
ground hard and essentially smooth during early summer. The boundary layer 
on the desert floor in early summer is thus quite similar to that on a smooth flat 
plate pO| . The measurements were made in the early summer season roughly 
between 6 PM and 9 PM during which nearly neutral stability conditions pre- 
vailed. Data set II was acquired over a rough terrain with ill-defined fetch at the 
meteorological tower at the Brookhaven National Laboratory. In sets I and II, 
data were acquired at heights 6 m and 35 m respectively. They were recorded 
simultaneously from two single hot-wire probes separated in the spanwise direc- 
tion y by 55 cm and 40 cm respectively. In both cases, the separation distance 
was within the inertial range, and was set nominally orthogonal to the mean 
wind direction (see below). Set III was acquired from an array of three cross- 
wires, arranged above each other at heights 11 cm, 27 cm and 54 cm respectively. 
These measurements are thus much closer to the desert floor. The hot-wires, 
about 0.7 mm in length and 6 /im in diameter, were calibrated just prior to 
mounting them on the meteorology towers and checked immediately after dis- 
mounting. The hot-wires were operated on DISA 55M01 constant-temperature 
anemometers. The frequency response of the hot-wires was typically good up to 
20 kHz. The voltages from the anemometers were suitably low-pass filtered and 
digitized. The voltages were constantly monitored on an oscilloscope to ensure 
that they did not exceed the digitizer limits. Also monitored on-line were power 
spectra from an HP 3561A Dynamic Signal Analyzer. The wind speed and di- 
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Figure 4: The one-dimensional energy spectrum computed from data set II. 
(From Ref. §].) 

rection were independently monitored by a direction indicator mounted on the 
tower (sets I and III), or a vane anemometer a few meters away (set II). The 
real-time durations of data records were limited only by the degree of constancy 
demanded of the wind speed and its direction. 

The Taylor microscale Reynolds number was 10,000 for set I and about 
20,000 for set II [||. For set III, the corresponding numbers were 900, 1400 
and 2100, respectively, at the three heights. Table ^ lists a few relevant facts 
about the data records analyzed here. Figure |^ shows a compensated third order 
longitudinal structure function and indicates the region over which the inertial 
range is expected to hold. As another example of the nature of the data, we 
show in Fig. ^ the longitudinal energy spectrum which displays extensive range 
of scaling. The slope is slightly larger than 5/3, as the compensated spectrum 
in the inset is supposed to clarify. 

4 Anisotropic contribution in the case of homo- 
geneity 

4.1 General remarks on the data 

We analyze data sets I and II described in section 3 (see Table |l|) to extract 
the lowest order anisotropic contribution to scaling (Q,Q). First, some pre- 
liminary tests and corrections need to be made. To test whether the separation 
between the two probes is indeed orthogonal to the mean wind, we computed 
the cross-correlation function {ui{t + T)u2{t)) . Here, ui and U2 refer to velocity 
fluctuations in the direction of the mean wind, for the arbitrarily numbered 
probes 1 and 2 respectively (see Fig. 0). If the separation were precisely orthog- 
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onal to the mean wind, the cross-correlation function should be maximum for 
r = 0. Instead, for set I, we found the maximum shifted shghtly to t = 0.022 s, 
implying that the separation was not precisely orthogonal to the mean wind. To 
correct for this effect, the data from the second probe were time-shifted by 0.022 
s. This amounts to a change in the actual value of the orthogonal distance. We 
computed this effective distance to be A « 54 cm (instead of the 55 cm that was 
set physically). For set II, the effective separation distance was estimated to be 
31 cm (instead of the physically set 40 cm). Next, we tested the isotropy of the 
flow for separations of the order of A. Define the "transverse" structure func- 
tion across A as S't(A) = {[ui{Ut) — U2{Ut)]'^) and the "longitudinal" structure 
function as S'l(A) = ([ui{Ut + UtA) ~ ux{Ut)\^) where iA = A/f7. If the flow 
were isotropic we would expect Q 

MA) = ^.(A) + y^§^. (30) 

In the isotropic state both longitudinal and transverse components scale with 
the same exponent, S't,l(A) oc A^^, and the ratio St/Sl is computed from (|o|) 
to be I-I-C2/2 ~ 1.35, because C2 ~ 0.69 (see below). The experimental ratio was 
found to be 1.32 for set I, indicating that the anisotropy at the scale A is small. 
This same ratio was about 1.8 for set III, indicating higher degree of anisotropy 
in that scale range. The differences between the two data sets seem attributable 
partly to differences in the terrain and other atmospheric conditions, and partly 
to the different distances from the ground. 

Lastly, we needed to assess the effects of high turbulence level on Taylor's 
hypothesis. A comparison was made of the structure functions of two signals 
with turbulence levels differing by a factor of 2 and no difference was found. The 
correction scheme proposed in |27| also showed no changes. For a few separation 
distances, the statistics of velocity increments from two probes separated along 
n, the direction of the mean wind, agreed with Taylor's hypothesis. More details 
can be found in Ref. §. 

The use of the isotropy equation ( |30| ) presents some evidence regarding the 
prevalence of modest amounts of anisotropy in the second-order statistics. This 
persistence is more explicit if one considers another familiar object, namely the 
cross-spectral density between horizontal and vertical velocity components. We 
consider the one-dimensional cross-spectral density (or shear-stress cospectrum) 
Ei^(k-i), which is zero in the case of isotropy. From dimensional considerations, 
the scaling exponent for this quantity is —7/3 (see [|l^ for theory and ||2^ for 
experimental tests). Figure ^ shows the cospectrum computed for the height 
of 0.54 m. The inset shows that the cospectrum compensated with a scaling 
exponent of —2.1 is flat. To the extent that this is numerically smaller than 
7/3, the decay of anisotropy is slower than expected for second-order quantities. 
Even allowing for the fact that the dimensional analysis assumes K41 scaling 
and therefore does not account for possible intermittency corrections in the 
anisotropic sectors, it is still clear that the small scales do not attain isotropy 
as fast as dimensional considerations suggest. We would like to use the proper- 
ties of the SO (3) decomposition in order to disentangle the isotropic from the 
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Figure 5: Log-log plot of the shear-stress cospectrum £13(^3) computed at 0.54 
m. The inset shows a log-log plot of kl'^Ei^^ks) vs. k^. The flat region indicates 
a region of scaling with exponent -2.1. 

anisotropic contributions and better quantify the anisotropy. We first consider 
the second-order structure function in detail. 

4.2 The tensor form for the second-order structure func- 
tion 

4.2.1 The anisotropic tensor component derived under the assump- 
tion of axisymmetry 

To obtain a theoretical form of the structure function tensor we first select a 
natural coordinate system. Our choice is to have the mean-wind direction n 
along the z-axis which we label as 3. The angle with respect to the 3-axis is 
the polar angle 9. The second axis is given by the separation vector A between 
the two probes. We make the simplifying assumption that the main symmetry 
broken in the flow is the cylindrical symmetry about the mean-wind direction. 
In other words, the main anisotropic contribution is cylindrically symmetric 
about the mean-wind direction. It is shown a posteriori that this assumption 
probably accounts for most of the anisotropy for this particular geometrical set- 
up. In the next section, we provide a complete analysis with no assumptions 
about the symmetries of the flow. We conclude from that analysis that the 
simplification of the present section is not inconsistent with the properties of 
the full tensor. 

Next, we write down the tensor form for the general second-order structure 
function (defined by Eq. ( ]l3| ) for n = 2) in terms of irreducible representations 
of the SO (3) rotation group. Since we are far from the wall, our interest is in 
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relatively modest anisotropies, and so we focus on the lowest order correction 
to the isotropic (j = 0) contribution. We therefore write 



S-^ir)^Sf£,{r)+Sfi,{r) + 



(31) 



We do not have a j = 1 term because it vanishes from parity considerations (the 
structure function itself is even in r in our experimental configuration) or by the 
incompressibility constraint. Assumed cylindrical symmetry of the anisotropic 
j = 2 contribution implies that we only consider the m = subspace in this 
contribution. Now, the most general form of the tensor can be written down by 
inspection. The case j = is well known (sec Eq. (pl|)) to be 



5"fo« = co(r) 



(2 + C2)<5"^-C2 



(32) 



where co(r) = cqt^^, and cq is a non-universal numerical coefhcient that needs 
to be obtained from fits to the data. The j — 2 component has six independent 
tensor forms and corresponding coefficients. These correspond to the different 
a^oo described in section 1. As with j = 0, the j = 2 component can be simplified 
by imposing conditions of incompressibility and orthogonality with the j = 
part of the tensor. This leaves us (in the case of cylindrical symmetry) with two 
independent coefficients which we call a and 6, giving, 



a(2) 



(Cf^ 



2)5' 



a/3 



C: 



2) 



2 -6) 
r"r^(n-r)2 



([Cf]^+3d^)+6)n"n^ 

id'' - 2 



^(2) 
S2 



(r"n^ + r'^n")(n-r) 



(2) 



(C 

(2Cf) 



-((2 

2) — 



(2) 



+ 3)(cf +2)<5"'^(n.r) 



a/3( 



(Cf^ 



-3)(cf^ +2)n"n'^ 



l)(Cf^ 



2) 



r"r^(n-r)" 



- ([Cf']^-4)(r-"n^ + r'^n")(n-r) 



(33) 



We note that the forms for the tensor were derived on the basis of the assump- 
tion of cylindrical symmetry. We used neither the Clebsch-Gordon method nor 
the "adding indices" method of section 1. This is why they were not automat- 
ically orthogonal to the j = contribution; that condition had to be explicitly 
imposed. However, in the end, this method also yields the same number of inde- 
pendent coefficients as the formal Clebsch-Gordon methods. This merely shows 
that there are different ways of writing the tensor forms for the basis elements. 



21 



However, the Clebsch- Gordon rules tell us the number, parity and symmetry of 
the terms we must in the end have for a given j, m sector. 

Finally, we note that in the present experimental set-up only the component 
of the velocity in the direction of n (the 3-axis) is measured. In the coordi- 
nate system chosen above we can read from Eqs. (^l|) and ( p3| ) the relevant 
component to be 



2 + C2 - C2 cos^ 9 



+ [(Cf +2)^-Cf^(3Cf +2)cos^^ 



.A 



(2) ^/■(2) 



2) cos^ 



(34) 



cos^ e 



(Cf +2)(cf +3)-Cf (3Cf +4) 
l)(Cf -2)cos4 0" . 



Here 9 is the angle between r and n, and r has been normalized by A, making all 
the coefficients dimensional, with units of (m/sec)^. Taylor's hypothesis allows 
us to obtain components of S"^ from (p^, with 9 = and variable, respectively. 
In other words, 



S^^{r, 9 = 0) = {[uiiUt + Utr) - ui{Ut)f) 
where = r/U, and 

533(r, 9) = {[ui{Ut + Utf) - U2{Ut)f). 



(35) 
(36) 



Here 9 = arctan(A/C/if ), tf = f/f/, and r = ,/A2 + {Uiff. 

The quantities on the left hand side of Eqs. (^ and (|3^) were computed 
from the experimental data and fitted to the theoretical expression ( p^ ) using 
the appropriate values of 9. The fits were performed in the range 1 < r/A < 10 
(0.54 m < r < 5.4 m) for set I and 1 < r/A < 25 (0.31 m < r < 8 m) for set II. 
The ranges were based on the constancy of the third order structure function 
(see Fig. ^ for an example). Panels (a) of Figs. || and Q show, for data sets I and 
II respectively, a comparison between the measured S^^ (r, 6* = 0) and the j = 
form of the equation. The comparison shows that the agreement is modest, and 
that the best-fit yields the exponent C2 to be 0.69. A careful analysis of the data 
elsewhere ||2^ shows that, if the effect of the shear is removed in a plausible way, 
the power-law fit is excellent over a range of scale separations. To include the 
j = 2 contribution, we fixed C,2 to be 0.69 and performed the following analysis. 

f2) 

For given values of the variables r and 9, we guessed the second exponent Q 
and estimated the unknown coefficients cq, a and h by using a linear regression 
algorithm. We followed this procedure repeatedly for different values of Q 
ranging from to 2. We then chose the value of (j^^ that minimized (the 
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C2 


/■(2) 
S2 


Co 


a 


b 


0.69 


1.38±0.15 


0.023±0.001 


-0.0051±0.0006 


0.0033±0.0005 


0.69 


1.36±0.10 


0.112±0.001 


-0.052±0.004 


0.050±0.004 



Table 2: The scaling exponents and the three coefficients in units of (m/sec)^ 
as determined from the nonlinear fit of Eq. (7) to data sets I (first line) and II 
(second fine). 

sum of the squares of the differences between the experimental data and the 
fitted values). 

In Fig. we present values as a function of (^"^ . The optimal value of this 

. . . (2) 

exponent and the uncertainty determined from these plots is Q ~ 1-38 ± 0.15 

(2) 

from set I, and Q ~ 1.36 ± 0.1 from set II. The best numerical values for the 
coefficients are presented in Table ||. Panels (b) in Figs. ^ and |7| show fits to 
the sum of the j = and j = 2 contributions to the experimental data. Even 
though the j = 2 contributions are small, they improve the fits tremendously. 
This situation lends support to the essential correctness of the present analysis. 
The figures show that the purely longitudinal structure function corresponding 
to 6* = is somewhat less affected by the anisotropy than is the finite 6 structure 
function (see especially, Fig. Q). The reason is the closeness of the numerical 
absolute values of the coefficients a and b (see Table |2|). For 6 = the two tensor 
forms multiplied by a and b coincide, and the j = 2 contribution becomes very 
small. The value of ^2 — 0.69 quoted above can be obtained from such a fit to 
the C2 = part alone; as long as one measures only this component, it seems 

reasonable to proceed with just that one exponent. However, the inclusion of 

(2) 

the second exponent Q improves the fit even for the longitudinal case; for the 
finite 9 case, this inclusion appears essential for a good fit. In fact, the fit using 
j = and j = 2 for data set II extends all the way to 40 m (beyond the range 
shown in Fig. |^). This indicates that the "inertial-range" where scaling theory 
applies is much longer than anticipated by traditional log- log plots. The close 
agreement with the theoretical expectation of 4/3 (e.g., Refs. ||l^,||l^), and the 
apparent reproducibility of the result for two different experiments is a strong 
indication that this exponent may be universal. 

It should be understood that, for the objects considered in this section, 
the exponent (^2^^ is just the smallest exponents in a possible hierarchy that 
characterizes higher order irreducible representations indexed by j. (The case of 
^2^^ will be discussed in section 5.) As discussed earlier, we expect the exponents 
to be a non-decreasing function of j, and that the highest values of j are being 
peeled off quickly when r decreases. Nevertheless, the lower order values of C^"''' 
can be measured and computed. The results of this section were first reported 
in 1^ . Results are also in agreement with the subsequent analysis of numerical 
simulations 
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(a) (b) 

Figure 6: The structure functions S'^^ for = and for non-zero 9 computed for 
set I. The dots are for experimental data and the hne is the analytic fit. Panel 
(a) presents fits to the j = component only, and panel (b) to components 
j = and j = 2 together. 
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Figure 7: The structure functions S^^ for = and for non-zero 9 computed 
for set III. The dots are for experimental data and the Une is the analytic fit. 
Panel (a) presents fits to the j = component only, and panel (b) to the sum 
of the components j = and j = 2. 
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Figure 8: The determination of the exponent ^2^'* from a least-square fit of 
S^^{r,9) to its analytic form. From set 1 we obtain a numerical value of the 
best fit, ^ = 1-38 ± 0.15 while from set III the best fit value is 1.36 ± 0.1 
both of which are in close agreement with the theoretical expectation of 4/3 of 
section 2.4 (without intermittency corrections). The differences in the nature 
of the minima are not understood. 



4.2.2 The complete j = 2 anisotropic contribution 

In the previous section, preliminary results on the scaling exponent (^2'^), were 
obtained under the assumption of cylindrical symmetry of the dominating anisotropic 
contribution. The analysis here is more complete, and takes into account the 
full tensorial structure. We show that this is feasible, and the final results are 
in agreement with those presented in the previous section. 

The method used to extract the unknown anisotropic scaling exponent is 
essentially the same as in the previous section. Since we look for the lowest order 
anisotropic contributions in our analyses, we perform a two-stage procedure to 
separate the various sectors. First we look at the small scale region of the inertial 
range to determine the extent of the fit with a single (isotropic) exponent. 
We then seek to extend this range by including appropriate anisotropic tensor 
contributions, and obtain the additional scaling exponents using least-squares 
fitting procedure. This procedure is self-consistent. 

Reference [|8| presents a detailed analysis of the consequences of Taylor's 
hypothesis on the basis of an exactly soluble model. It also proposes ways for 
minimizing the systematic errors introduced by the use of Taylor's hypothesis. 
In light of that analysis we will use an "effective" wind, UcS, for surrogating 
the time data. This velocity is a combination of the mean wind U and the 
root- mean-square u', 

L/cff = \/u^ + {hu'Y (37) 
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where 6 is a dimensionless parameter. 

In the second order structure function defined already, viz., 

S'"^(r) = ((m"(x + r) - u"(x))(w''(x + r) - ^'^(x))), 



(38) 



the j = 2 component of the SO (3) symmetry group corresponds to the lowest 
order anisotropic contribution that is symmetric in the indices, and has even 
parity in r (due to homogeneity). Although the assumption of axisymmetry 
used in Q and in the previous section seemed to be justified from the excellent 
qualities of fits obtained, we attempt to fit the same data (set I) with the 
full tensor form for the j = 2 contribution. The derivation of the full ) = 
2 contribution to the symmetric, even parity, structure function appears in 
appendix A. We then find the range of scales over which the structure function 



S'^'-'{r,e = 0) = {{u[^\x + r) 



,(3) 



(39) 



with the subscript 1 denoting one of the two probes, can be fitted with a single 
exponent. To find the j — 2 anisotropic exponent we need to use data taken 
from both probes. To clarify the procedure, for the geometry shown in Fig. ||, 
what is computed is actually 



S'\r, 9) = {[uf\u,iit + C/efft,) - uf{U,st)f). 



(40) 



Here 9 — arctan(A/C/offif ), tf — r/Ucs, and r = + (Ocffif )^- C^cff is 

defined by Eq. ( p7[ ) with the optimal value of b taken from model studies to be 
3. For simplicity we shall refer from now on to such quantities as 



533(r,0) = {{uf'>{x + r)-u'i\x))^) 



(41) 



Next, we fix the scaling exponent of the isotropic sector to be 0.69 and find 
the j = 2 anisotropic exponent that results from fitting to the full j — 2 tensor 
contribution. Wc fit the objects in Eqs. ( |39| ) and ( |4l| ) to the sum of the j = 
(with scaling exponent ^2 = 0-69) and the j = 2 contributions. The sum is 
given by (see appendix A) 



S'33(r, 9) 



S'ioir,0)+Sf^,{r,9) 
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C2 r 



2 + C2 - C2 cos^ 
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(2) 



(d^)+2)^-Cr^(3C^^^+2)cos^ 



(2)/qA(2) 
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Figure 9: The structure functions computed from data set I and fit with the 
j = and full j = 2 tensor contributions using the best fit values of exponents 

(2) 

C,2 = 0.68 and Q = 1-38. Panel (a) shows the fit to the single-probe (9 = 0) 
structure function in the range 0.2 < r/A < 25 and panel (b) shows the fit to 
the ^-dependent structure function in the range 1 < r/A < 25. 



^9,2,2 



+ ai,2,2 



■aJ 



+ 2C:^^'(Cf' - 2)cos3 6lsin( 
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2Cf(Cf)-2)sin^ 



We fit the experimentally generated functions to the above form using values 

(2) 

of Q ranging from 0.5 to 3. Each iteration of the fitting procedure involves 

f2) 

solving for the six unknown, non-universal coefficients. The best value of Q 
is the one that minimizes the value for these fits; we obtain that to be 
1.38 ± 0.15. The fits with this choice of exponent are displayed in Fig. ||. 

The corresponding values of the six fitted cocfhcients is given in Table 2. The 
range of scales that are fitted to this expression is 0.2 < r/A < 25 for the 6' = 
(single-probe) structure function and 1 < r/A < 25 for the 6* 7^ (two- probe) 
structure function. We are unable to fit Eq. ( ^ ) to scales larger than about 
12 meters without losing the quality of the fit in the small scales. This limit 
is roughly twice the height of the probe from the ground. Based on the earlier 
discussion, we should be in the regime of the largest scales up to which the 
three-dimensional theory would hold. (Beyond this scale the boundary layer is 
dominated by the sloshing motions which are quasi two-dimensional in nature.) 
Therefore this limit to the fitting range is consistent with our expectations for 
the maximum scale of three-dimensional turbulence. We conclude that the 
structure functions exhibit scaling behavior over the whole scaling range, but 
this important fact is missed if one does not consider a superposition of the 
2 = and j = 2 contributions. We thus conclude that the estimate for the 



2 scaling exponent C. 



(2) 



1.38. This same estimate was obtained in 
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and in the previous section using only the axisymmetric terms. The value of 
the coefficients a and h are again close in magnitude but opposite in sign — ^just 
as in [Q, giving a small contribution to S^^{r,9 = 0). The non-axisymmetric 
contributions vanish in the case of 6* = 0. The contribution of these terms to 
the finite 9 function is relatively small because the angular dependence appears 
as sin 9 and sm^9, both of which are small for small 9 (large r); it was for 
this reason that we were able to obtain in the previous section a good fit to 
just the axisymmetric contribution. Lastly, we note that the total number of 
free parameters in this fit is 7 (6 coefficients and 1 exponent). The relative 
"flatness" of the function near its minimum (see, especially the left panel 
of Fig. 1^) may be indicative of the large number of free parameters in the fit. 
However, the value of the exponent is perfectly in agreement with the analysis 
of numerical simulations j^] in which one can properly integrate the structure 
function against the basis functions, eliminating all contributions except that of 

(2) 

the j = 2 sector. Furthermore, fits to the data in the vicinity of Q = 1.38 show 
enough divergence from experiment that we are satisfied about the genuineness 
of the result. The results of this analysis taking into account the full j = 2 
tensor were first presented in [l^ . 

4.3 Summary 

We consider the lowest order anisotropic contribution to the second order tensor 
functions of velocity in the atmospheric boundary layer. The j = 1 contribu- 
tion is absent in our experimental configuration because of the incompressibility 
constraint. The j ~ 2 sector is expected to be the dominant contribution to 
anisotropy. First we make the a priori assumption that the cylindrical symmetry 
is broken in the first deviation from anisotropy. The use of the Clebsch-Gordon 
rules tell us the number of terms that must be included. We derive on this basis 
the tensor form for the j = 2 contribution that is axisymmetric. The conditions 
of orthogonality with j = and incompressibility are used to constrain the un- 
known coefficients giving two unknown coefficients and one unknown exponent 
in the anisotropic sector. Using this form along with the known isotropic contri- 
bution, we can extract the anisotropic scaling exponent from the experimental 
data. We have also used the full form of the j = 2 tensor including all its 2j + 1 
components in order to examine if the results from the initial analysis were jus- 
tified. The excellent agreement between the two strengthens our confidence in 

(2) 

the value of the j = 2 exponent Q ~ 1.36, and in our initial assumption of the 
breaking of axisymmetry of the flow. 

5 Anisotropic contribution in the case of inho- 
mogeneity 
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5.1 Extracting the j=l component 

The homogeneous structure function defined in Eq. ( ^ ) is known from proper- 
ties of symmetry and parity to possess no contribution from the j — 1 sector 
(see section B.2), the j — 2 sector being its lowest order anisotropic contributor. 
In order to isolate the scaling behavior of the j = 1 contribution in atmospheric 
shear flows we must either (a) explicitly construct a new tensor object which will 
allow for such a contribution, or (6) extract it from the structure function itself 
computed in the presence of inhomogeneity. Adopting the former approach, we 
construct the tensor 

r"'3(r) = (u"(x + r) - u"{x)){u^{yL + r) + u'^(x))). (43) 

This object vanishes both when a — /3 and when r is in the direction of ho- 
mogeneity, viz., the strcamwise direction. From data set III we can calculate 
this function for non-homogeneous scale-separations (in the shear direction) . In 
general, this will exhibit mixed parity and symmetry; the incompressibility con- 
dition does not reduce our parameter space. Therefore, to minimize the final 
number of fitting parameters, we examine only the antisymmetric contribution. 
In section B.2, we have derived the antisymmetric tensor contributions in the 
j = 1 sector, and used this to fit for the unknown j = 1 exponent. We describe 
the results of this effort below. This can be used to find j = 1 exponent for the 
inhomogeneous structure function which is symmetric but has mixed parity. We 
do not present the result of that analysis here essentially because it is consistent 
with those from the antisymmetric case. 

Returning to consideration of the antisymmetric part of the tensor object 
defined in Eq. (p3|), viz., 

^ r"^(r)-r/^"(r) 

^ ' 2 

= (u"(x)u''(x + r)) - (u'^(x)u"(x + r)), (44) 

it is easy to see that it will only have contributions from the antisymmetric 
j = 1 basis tensors. An additional useful property of this object is that, for the 
configuration of data set II, it does not have any contribution from the isotropic 
J = sector spanned by 5°^^ and r°'r^ since these objects are symmetric in 
the indices. This allows us to isolate the j = 1 contribution and determine its 
scaling exponent ^2^^ starting from the smallest scales available. Using data (set 
III) from the probes at 0.27 m (probe 1) and at 0.11 m (probe 2) we calculate 

r3i(r) = {ui\ic)ui\i^ + v)) - (wf^(x + r)4'^(x)) (45) 

where again superscripts denote the velocity component and subscripts denote 
the probe by which this component is measured. The goal is to fit this experi- 
mental object to the tensor form derived in appendix B, Eq. (^7|), namely, 

f^^{r,e,(t) = 0) = -03,1^0'^''^" sin 6* + a2,i^ir'^^^'^ + 034 -ir'^^'' cos 6*. (46) 
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Figure 10: The minimization by the best-fit value of the exponent Cj^^ of the 
j = I anisotropic sector from the fit to ^-dependent T^^ {r, 9) function in the 
range 1 < r/A < 2.2. 







a2,i,i 


as. 1,-1 


1 ±0.15 


0.0116 ±0.001 


0.0124 ±0.001 


-0.0062 ±0.001 



Table 3: The values of the exponents and coefficients (in units of (m/sec)^) 
obtained from the fit to the function T'^^(r, 0). 



Figure |10| gives the x minimization of the fit as a function of Q 



(1) 



We 



obtain the best value to be 1 ± 0.15 for the final fit. This is shown in Fig. 

The fit in Fig. 11 peels off at around r/A = 2. The values of the coefficients 
corresponding to the exponent = 1 a-re given in Table ^ The maximum 
range of scales over which the fit works is of the order of the height of the 
probes from the ground, consistent with the considerations presented earlier. 
This value of the scaling exponent of the j = 1 sector is in agreement with 
the theoretically expected value of unity (section 2.4). Again, we have satisfied 
ourselves that a different value of the exponent yields a substantially poorer fit 
to the data. These findings significantly strengthen the proposition |^ that the 
scaling exponents in the various sectors (at least up to j = 2) are likely to be 
universal. 



6 The higher order structure functions 
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Figure 11: The fitted T'^^(r, 6') function. Tlie dots indicate tlie data and tlie 
line is the fit. 



6.1 Introduction 

The method of SO (3) decomposition is applicable to tensor functions of any 
rank. However, computing the anisotropic terms in the higher order structure 
functions (n > 2) becomes an increasingly difficult task. The number of con- 
tributions to a given j (i.e., the q index) quickly becomes too large. If we were 
to use the methods adopted in the case of the second order structure function, 
the large number of free parameters available for fitting the data guarantees 
over-fitting. This would not allow us to unambiguously extract the anisotropic 
scaling behavior. In order to circumvent this problem we recall two proper- 
ties of the tensor decomposition. First, for the single-point measurements, with 
= (no angular dependence), although the tensorial information is lost, the 



free coefficients collapse to a single number, see Eqs. ( |29| ) and (28). Second, 
the isotropic part of the structure function tensor is very easily obtained, as 
we will show below. Examination of this part of the S0(3) expansion will show 
us which tensor components do not contribute to the isotropic sector. We then 
extract anisotropic exponents by considering only those tensor components that 
are explicitly zero in the isotropic sector, so that whatever is measured derives 
its contribution entirely from the anisotropic sector. We use an interpolation 
formula to compensate for the large-scale encroachment of inertial-range scales. 
This allows us to examine the lowest order anisotropic scaling behavior. The 
resulting anisotropic exponents for a given tensorial order are larger than those 
known for the corresponding isotropic part. One conclusion that emerges is that 
the anisotropy effects diminish with decreasing scale, although more slowly than 
previously thought. 

We can use the present method in principle to examine the anisotropic contri- 
bution of tensors of any order without requiring the knowledge of the particular 
mathematical form of the anisotropic sectors of these tensors. This is a consid- 
erable advantage theoretically because the high-order tensors are non-trivial to 
compute; it is an advantage experimentally because, unlike in numerical simula- 
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tions, one can measure only some components for simple geometric arrangements 
of probes. The results of this analysis were first presented in |13{. 



6.2 Method and results 

In this part of the analysis we use data set II described in section 2. We first 
consider the second-order tensor ^"^(r). Isotropy implies that this tensor can be 
expressed as a linear combination of two terms, 6"^ and r°'r^ , As is well known, 
both terms give non-zero contributions to longitudinal as well as transverse 
components, corresponding to a = (3. For a ^ P these two terms are identically 
zero if r is taken to be in the streamwise direction 3. Therefore, we compute 
the so-called mixed structure function 

53i(r) = + r) - u'^{x)){u^{x + r) - u^{x))), (47) 

where, as already noted, the superscripts 1 and 3 denote the vertical and stream- 
wise components respectively. This object is identically zero in the isotropic 
sector, and so, any non-zero value comes from anisotropy. In particular, any 
scaling behavior that it obeys should relate solely to anisotropy. By computing 
Eq. ( ^7| ) and examining its scaling, we intend to extract the purely anisotropic 
scaling behavior in the j = 2 sector, uncontaminated by any isotropic scaling, 
in contrast to the case of either longitudinal or transverse structure functions. 



6.2.1 The second-order structure function 

The previous paragraph provides the motivation for examining the measured 
structure functions S^^{r). However, as we shall see shortly, apart from the 
expected behavior in the dissipative range and saturation at some large scale, 
there appears to be no distinct inertial range scaling. We suspect that this 
happens because there is poor scale separation, since the probes are fairly close 
to the ground; in fact, the large scales (which we expect to be of the order of the 
height of the probe from the ground and larger (see previous section)) may be 
encroaching significantly into the inertial range. We would be aided materially 
in our search for scaling if, somehow, the large-scale effects can be separated. 
One way of doing this is to write down an interpolation function that models 
the entire structure function in its three different scaling regions — a dissipative 
range that scales like when r is of the order of the Kolmogorov scale rj, a 
large-scale behavior that tends to saturate (indicating decorrelation) as r gets 
to be larger than L, and the intermediate inertial range for rj « r « L which 
may exhibit scaling. Through the use of the interpolation formula, one can 
extract the scaling part in a natural way. This is described below. 

A suitable form of the interpolation function is given in Q for structure 
functions of arbitrary order. It has the form 
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r(m) 



Figure 12: Log-log plots of transverse structure functions at 0.54m. x denote the 
second order, O the fourth order and the solid lines represent the interpolation 
fit. 



where A„, _B„, C„ and D„ are variable parameters. This formula is an extension 
of that given in Ref. psf and includes a large-scale term. Such extensions have 
been attempted earlier (e.g., Ref. jljl). Antonia et al. Q have successfully tested 
the function (psf). Dhruva |^ has shown that the present interpolation formula 
works extremely well for longitudinal structure functions of order 2, 4 and 6. To 
reinforce this point, we test its performance by comparing it to the measured 
transverse structure function, a — f] = 1, r in the direction 3. For each data set, 
the height of the probe is assumed to be the large scale L. The fit is shown for 
the transverse structure function of orders 2 and 4 at the 0.54 m probe in Fig. 
The agreement between the formula and the data is excellent. Taken together 
with similar conclusions in |6j for longitudinal structure functions, we conclude 
that the interpolation formula describes the familiar structure functions very 
well. For this pragmatic reason, we shall adopt it for our purposes here, and 
test the robustness of the results obtained in the appendix C. 

In the formula (^8|), the large scale behavior is given by the factor (1 -|- 
D2{r/L))'^^^^^. If the measured structure function is divided by this factor, 
we should recover the contribution of the remaining parts — in particular the 
inertial range part, with the leading order scaling exponent given by 2 — 2C2. 
Figures |l^ and |lj display a second-order anisotropic structure function for two 
heights above the ground. Presumably because of the finiteness of the Reynolds 
number and the relatively large shear effect, the scaling in the intermediate 
range ry << r << L is not apparent. However, by dividing out the large 
scale contribution as described above, we see two distinct regions of scaling; the 
dissipative range of and the extended mid-range which scales with exponent 
between 1.22 and 1.12. The advantage of the scaHng function is thus evident: 
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Figure 13: Log-log plot of second order mixed structure function at 0.54m. X 
denote data, the solid line is the interpolation fit (not visible beyond an r of 
10~"'^m because of the closely packed symbols), and O correspond to the large- 
scale compensated function. 
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it has allowed us to extract a scaling exponent that is most likely to be due to 
anisotropy. The values of the fitted parameters and the corresponding Q are 
given in Tables ^ and || for the probes at 0.54 m and at 0.27 m respectively. The 

error on the measurement of the exponent C2 at 0.54 m is about 0.05 while at 

(2) 

0.27 m it is about 0.08. This gives an error on the estimates of Q of O-O''' ^-nd 
0.11 respectively, and places the theoretically expected value of « 4/3 within 
1.5 to 2 standard deviations of the present value. This result is consistent with 



general expectations 1 15 and the findings of the previous section. 



6.2.2 Higher-order structure functions 

In general, the tensor forms contributing to the j = sector for tensors of 
any rank n are composed of linear combinations of the Kronecker-(5 and the 
components of r along the tensor indices. The following is a list of isotropic 
tensor contributions for rank 3 through 6: 

• n — 3: S"^r^ + permutations, and r'^r^r^; 

• n — 4: 5°'^5'^^ + permutations, S"^r''r^ + permutations, and 

• n — 5: S"^S''^r'^ + permutations, S'^^r^r^r'^ + permutations, and r^r^r^r^r^; 

• 71 = 6: S^^S^^S^'^ + permutations, S^'^S^^r'^r'^ + permutations, 
gaPj.f^S^ti^i' ^ permutations, and r°'r^r''r^r'^r'^ . 

Based on the above considerations, it can be expected that the structure 
function components that are zero in the j — sector are: 

• 71 = 3: 5"^ (transverse), S^^^; 

.71 = 4: ^3331^^3111. 

.71 = 5: (transverse), 533111^ 533331. 

. 71 = 6- 5'333111 ^311111 ^333331 

Note that the odd-order transverse structure function is always zero in the 
isotropic sector. The functions we shall now consider are given in the second 
column of Tables || and ||. For the case of the third and fifth order transverse 
structure functions we use the moments of the absolute value of the velocity 
differences in order to obtain better convergence. In using the interpolation 
function we assume that the inertial range scaling of these anisotropic com- 
ponents is given by a single exponent (n^ where the superscript denotes an 
isotropic exponent without reference to the precise j sector. The compensated 

functions (with large-scale effects removed) are shown in Figs. |l5|-^. The errors 

(2) 

on the value of Qi obtained are about 7% at 0.54 m and about 9% at 0.27 m. 
For comparison, the last column in Tables ^ and || gives the empirical values of 
the isotropic scaling exponent of the same order as given in The entries in 
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Order n 


Tensor 




Bn 




Dn 




Cn 


2 




3.9 


0.014 


0.39 


0.67 


1.22 


0.7 


3 




2400 


0.010 


0.93 


2.28 


1.14 


1 


4 


^3331 


5200 


0.014 


1.21 


0.27 


1.58 


1.26 


5 


^11111 


1.22 X 10^ 


0.029 


1.59 


3.09 


1.82 


1.56 


6 


^333111 


3.75 X 10^ 


0.041 


1.93 


0.50 


2.14 


1.71 



Table 4: Structure function calculated and the anisotropic scaling exponents for 
the data at 0.54m. 



Order n 


Tensor 


An 


Bn 


Cn 


Dn 


Cn — 2Cn 


Cn 


2 


^31 


9.4 


0.005 


0.44 


0.52 


1.12 


0.7 


3 




6940 


0.015 


0.89 


2.78 


1.21 


1 


4 


^3331 


2.1 X 10^ 


0.014 


1.23 


0.23 


1.54 


1.26 


5 


^11111 


5.9 X 10^ 


0.028 


1.58 


3.52 


1.84 


1.56 


6 


^333111 


2.7 X 10^ 


0.038 


2.00 


0.34 


2.00 


1.71 



Table 5: Structure function calculated and the anisotropic scaling exponents for 
the data at 0.27m. 



this column are measurably smaller than the corresponding non-isotropic expo- 
nents. This suggests that the isotropic component alone survives at very small 
scales. 

6.3 Summary 

We have presented a new method of extracting anisotropic exponents that avoids 
mixing with contributions from the isotropic sector. We do this by explicitly 
constructing those tensors that are zero in the isotropic sector. An operational 
step in the extraction of the scaling exponents is the use of an interpolation 
formula in the spirit of a "scaling function". This method has allowed us to 
examine anisotropic effects in structure function tensors of order greater than 
2. The resulting anisotropic exponents are consistently larger than those known 
for isotropic parts at all orders. This strongly suggests that anisotropy effects 
decrease with decreasing scale. However, the rate of decrease is much slower 
than expected from dimensional arguments (which yield 4/3, 5/3, 2, 7/3 and 
8/3 for orders 2 through 6, to be compared with the values obtained at 0.54 m 
of 1.22, 1.14, 1.58, 1.82, 2.14). 

Our conclusions are based on the use of the interpolation formula, given by 
Eq. (^). While this formula is not based on a solid theoretical framework, 
we have shown that it works very well in describing the measured structure 
functions. We have also performed tests of its robustness by fitting it to smaller 
sections of the data in order to detect changes in the exponent. A discussion 
of these checks and their results is presented in appendix C. To the lowest 
order, the results are independent of the r-segment to which the formula is 
fitted (except, of course, when the fit is entirely for the dissipation range or 
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Figure 15: Log-log plot of third-order transverse structure function at 0.54ni. 
The legend is the same as for Fig. |l^. 
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Figure 20: Log-log plot of fifth-order transverse structure function at 0.27m. 
The legend is the same as for Fig. 
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in the large-scale range). Any other formula that works equally well will yield 
similar results. Even so, the formula is empirical, which is why we have not 
paid much attention to the fact that the scaling exponents obtained for the two 
probe positions are slightly different, and that the second-order exponent for 
0.54 m is slightly larger than that obtained for the third-order. On the whole, 
the trend is that the anisotropic exponents become larger for larger orders of 
the structure function. 

An interesting conclusion of the present work is that the effects of anisotropy 
vanish with decreasing scale more slowly than expected. The expectation in the 
light of the SO (3) formalism is that a hierarchy of increasingly larger exponents, 
corresponding to increasingly higher-order anisotropic sectors 0, would exist. 
This expectation appears to be true in the case of the passively advected vector 
field where a discrete spectrum of anisotropic scaling exponents is obtained 
theoretically for all anisotropic sectors. In the present experiments, the fact that 
anisotropic effects can be fitted reasonably well by power laws suggests that the 
high-order effects may be small. It is perhaps true, however, that the power- 
laws described here may contain high-order corrections, and that the exponents 
deduced for the behavior of anisotropy may indeed undergo some revision when 
contributions from other sectors of the S0(3) decomposition are also consid- 
ered. In spite of this possibility, we wish to emphasize that the anisotropy 
effects for each order of the structure function appear to be well described by 
something close to a power law with a single exponent. The magnitudes of the 
anisotropic exponents in each order indicate that the roll-off from isotropy hap- 
pens less sharply than previously thought, but the roll-off occurs nevertheless. 
The higher-order objects considered here have not been studied extensively in 
the light of anisotropy. 

7 Conclusions 

We have used the S0(3) decomposition of tensorial objects such as structure 
functions measured in fluid turbulence experiments. This enables us to separate, 
for any given order structure function, the isotropic part from the anisotropic 
part. We have used experimental data at very high Reynolds numbers (Taylor 
microscale Reynolds numbers of up to 20,000) to carry through this decompo- 
sition. 

We considered the lowest order anisotropic contribution to the second order 
tensor functions of velocity in the atmospheric boundary layer. For second order 
structure functions, we have extracted the anisotropic scaling exponent from the 
experimental data. We also used the full form of the j = 2 tensor, as well as the 
axisymmetric version of it. In both instances, we have shown that the j = 2 ex- 

(2) 

ponent Q ' w 1.36. This compares favorably with the estimate based on purely 
dimensional grounds. It should be noted that the dimensional estimates do not 
take into account possible anomalous scaling. We expect that turbulence flows 
will exhibit anomalous scaling in the isotropic sector, and in every anisotropic 
sector in the SO (3) hierarchy. By suitably arranging the measurement configu- 
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ration, we have also extracted the j = 1 contribution, and find the appropriate 
scahng exponent to be close to the dimensional estimate of unity. 

For high orders, we have presented a new method for the extraction of 
anisotropic exponents that avoids mixing with contributions from the isotropic 
sector. We do this by explicitly constructing those tensors that are zero in 
the isotropic sector. The resulting anisotropic exponents are consistently larger 
than those known for isotropic parts at all orders. This strongly suggests that 
anisotropy effects decrease with decreasing scale. However, the rate of decrease 
is slower than expected from dimensional arguments (which yield 4/3, 5/3, 2, 
7/3 and 8/3 for orders 2 through 6, to be compared with the values obtained at 
0.54 m of 1.22, 1.14, 1.58, 1.82, 2.14). 

Armed with these details, it is now possible to say something of value about 
local isotropy. These statements are interspersed through out the text, but the 
most important conclusion is that the effects of anisotropy within a structure 
function of a given order do vanish with decreasing scale, though more slowly 
than expected. This provides a rich perspective on the notion of local isotropy 
and removes a hurdle towards the search for a universal theory of small-scale 
turbulence. 



A Full form for the j 
homogeneous case 



2 contribution for the 



Each index j in the SO (3) decomposition of an rt-rank tensor labels a 2j + 1 
dimensional S0(3) representation. Each dimension is labeled by m = —j, ~j + 
The j = sector is the isotropic contribution while higher order j's 
should describe any anisotropy. The j — terms are well known to be 



{2 + C2)S'-^ 



(49) 



where ^2 ~ 0.69 is the known universal scaling exponent for the isotropic con- 
tribution, and Co is an unknown coefficient that depends on the boundary con- 
ditions of the fiow. For the j — 2 sector, which is the lowest contribution to 
anisotropy to the homogeneous structure function, the m — (axisymmetric) 
terms were derived from constraints of symmetry, even parity (because of ho- 
mogeneity) and incompressibility on the second order structure function to be 



+ 2Cf (Cf -2) 



(Cf^-2)<5"'3-Cf (Cf^+6) 
r"r'^(n • r)^ 



+ ([Cf¥ + 3Cf +6)n°n'' 



(50) 
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(2) 



(Cf^+3)(cf^+2)5"^(n.rf 



a f) 



(Cf^ + 3)(Cf ' + 2)n"n'5 + (2Cf ^ + - 2) 

- i[C2^f - 4)(r"n'' + r^n")(n • r) 

(2) 

Here Q is the universal scaling exponent for the j — 2 anisotropic sector and a 
and b are independent unknown coefRcients to be determined by the boundary 
conditions. We would now like to derive the remaining m = ±1, and m = ±2 
components 

= E««.2,™r?^''B,l,„(f), (51) 
q 

where C| is the scaling exponent of the j = 2 S0(3) representation of the n = 2 
rank correlation function. The B'^j^^{r) are the basis functions in the S0(3) 
representation of the structure function. The q label denotes the different possi- 
ble ways of arriving at the the same j and runs over all such terms with the same 
parity and symmetry (a consequence of homogeneity and hence the constraint 
of incompressibility) ; see Ref. In our case, even parity and symmetric in 
the two indices. In all that follows, we work closely with the procedure outlined 
in [||. Following the convention in [|| the q's to sum over are q = {1,7,9,5}. 
The incompressibility condition daU°' = coupled with homogeneity can be 
used to give relations between the CLqj^m for a given {j,m). That is, for j — 2, 
m = —2 ... 2, we have 

(Cf ^ - 2)ai,2,™ + 2(cf ^ - 2)ar,2,m + (Cf ^ + 2)09,2,™ = 

ai.2,m + (d^^ + 3)07,2, m + Ci^^a5.2,m = 0. (52) 

We solve equations (3) in order to obtain 05, 2, m and 07,2,™ in terms of linear 
combinations of 

ai,2,,n([Cf ]^ - Cf' - 2) + a9.2,m{Ki'^? + 5cf ^ + 6) 
2Cf(Cf-2) 

ai,2,m(2 - C2^'') - a9,2,m(2 + C2^^) 
a7.2,m = TIT) ■ (03j 

2(cf^-2) 

Using the above constraints on the coefRcients, we are now left with a linear 
combination of just two linearly independent tensor forms for each m: 

S-t2,^ = «9,2,™rC^''[-Cf (2 + CfV7l™(f) 
+ 2Cf)(Cf -2)i?9t„(f) 

+ ([Cf¥ + 5Cf +6)S5i,™(f)] 



^,2 J -r us,2 -r ^)J^^a,ra 

ai,2,™rC^''[2Cf^(Cf -2)<,„(f) 
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-Cf)(Cf^-2)B«j,™(f) (54) 

+ ([Cf]'-Cf -2)B,t™(f)]. 

The task remains to find the explicit form of the basis tensor functions 



^",2.m(f), 9 e {1,7,9,5}, TO G {±1,±2} 



-'9,2,771 



We obtain the m = {±1, ±2} basis functions in the following derivation. We 
first note that it is more convenient to form a real basis from the r^l2?n(f) since 
we ultimately wish to fit to real quantities and extract real best-fit parameters. 
We therefore form the r^y2fe(f) (A; = —1, 0, 1) as follows: 



r'?2-i(f) 



r2y2o(r) =r^cos2 6i = r3^ 
oy2-i(r)-y2+i(?) 



r 

T 



(cos ^ — i sin 0) cos 9 sin 
-|-(cos ^ + i sin (f)) cos 6 sin 9 



r cos 9 sm U cos = r^ri 
,y2-i(f) + F2+i(r) 



-2i 



-2i 

(cos 

— (cos <j) + isin cj)) cos 9 sin 9 



isiuffl cost/smi 



= r cos0sin0sin(/) = r3r2 



r'F2-2(r) = 



,y22(f)-y2-2(r) 

2i 



r 
2i 



(cos 2(j) + i sin 2<f>) sin^ — (cos 2(j) — i sin 2(/)) sin^ 9 



r'l2+2(r) 



r'^ sin 2(/)sin^ 9 = 2rir2 

2y22(r) + y2-2(f) 



r 
~2 



(cos 20 + i sin 20) sin^ 9 + (cos 20 — i sin 20) sin^ 



= cos20sin^ ^ = rf — rf . 



(55) 



This new basis of r^y2fc(r) is equivalent to the r^Yjm{j^) themselves as they 
form a complete, orthogonal set (in the new k's). We omit the normalization 
constants for the spherical harmonics for notational convenience. The subscripts 
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on r denote its components along the 1 (to), 2 (p) and 3 (n) directions; m 
denotes the shear direction, p the horizontal direction parallel to the boundary 
and orthogonal to the mean wind direction and n the direction of the mean 
wind. This notation simplifies the derivatives when we form the different basis 
tensors; the only thing to remember is that 



= a"(r • m) = m" 
9"r2 = 9"(r-p) = p" 
a"r3 = (9"(r • n) = n". 



(56) 



We use the above identities to proceed to derive the basis tensor functions 
= r-2^"''(r .n)(r-m) 

= r-2[(r"m'3 + r^^m°'){r ■ n) + (r^n'^ + rl^n"){v ■ m)] 



B?iAr 

Da/3 (~ 



Q/3 

1,2,-2 



Q/3 



7,2,-2V^ 

a/3 Cf 
9,2,-2^,^ 

^5,l-2(f 

<,2(f 

57l2(f 



9,2,2 



Da/3 /- 
"5,2,2U 



'^(r • n)(r • m) 



r-25"''(r .n)(r-p) 

r-2[(r°/ + rV)(r ■ n) + {r°'n>^ + r'^n")(r ■ p)] 



_ j,-2^a^/3 



(r •n)(r-p) 



n"/ + n V 
2r-2^"'3(r .m)(r-p) 

2r-2[(rV + rV)(r ' m) + {r"m'^ + r^m°'){r ■ p)] 
2r"V"r'^(r-m)(r-p) 

2(mV + "^V) 
r-25"''[(r.m)2-(r-p)2] 

2r-2[(r-"TO'3 + r^TO")(r • m) - (rV + r-V)(r • p)] 



_ j,-2^a^/3 



[(r.m)2-(r.p)2] 



(57) 



Substituting these tensors forms into Eq. [5J we obtain the full tensor forms 
for the j = 2 non-axisymmetric terms, with two independent coefficients for 
each k: 



([C2 T + 5Cr' + 6)(n°TO'^ + n^m") 
+ 2Cf '(Cf ^ - 2)r-V"r'3(r • n)(r • m) 
- cf ^(2 + Ci^^)r-^[{r"m^ + r^m"){r ■ n) 

+ (r"n'' + r'3n") (r-m)] 
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+ ai,2,-ir^ 



"5^=2, fc=2 W 



+ 2Cf '(Cr^ - 2)r-2^«'5(r . n)(r • m) 

- Cf ^ (Cf ^ - 2)r-2[(r"m'^ + r^m«)(r • n) 



+ (r"n^ + r'^n")(r-m)] 

^([Cf^]' + + 6)(nV + nV) 



2)r 



^(r • n)(r • p) 



- (2 + Cf + A")(r • n) 

+ (r"n'' + r''^n")(r-p)] 

+ ai.a.ir'^^" [([Cf ^ " 2)(nV + nV) 
+ 2Cf (Cf -2K25"/3(r.n)(r-p) 

- Cf^(d'^ - 2)r-2[(rV + rV)(r • n) 



A2) 



{[(2? + ^2 + 6)(mV + mV) 



+ 2Cf '(Cf ^ - 2)r-Vr'^(r • p)(r • m) 
- 2Cf '(2 + V"^[(^V + rV)(r • m) 

+ (r"m'^+r'^m")(r-p)] 



a(2) 



^2([Cf ]'-Cf -2)(mV + mV) 
+ 2Cf (Cf -2)r-25"^'^(r.m)(r.p) 

- 2C^^) (Cf ' - 2)r-2 [(r"/ + r V ) (r • m) 

+ (r"m'3 + /^m") (r-p)] 

= a9,2,2r'^^" [2([Cf ¥ + 5Cf ^ + 6)(m"m^' - pV) 
+ 2Cf )(Cf ' - 2)r-V"r''^[(r • m)^ - (r • p)^] 

- 2Cf ^(2 + C^^^)r-^[{r"m'^ + r^m")(r • m) 

- (rV + rV)(r-p)] 



+ ai,2,2r' 



[2([Cf]^-d^^-2)(m«m'5-pV) 

+ 2Cf ^ (Cf ^ - 2)r-2^«'5 [(r • m)^ - (r • p)^] 
- 2Cf '(Cf ^ - 2)r-2[(r«m^ + r^m''){r ■ m) 
- (rV + rV)(r -P)] • 



(58) 
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4> = 7r/2,Q = /3 = 3 


= O,Q = ^ = 3 


= 0,a = ;3 = 1 


= 0,a = 3,;9 = 1 


k 




= 


5^ 


61 = 


e^o 


f = 


e 


61 = 





2 


2 


2 


2 


2 


2 


2 





-1 








1 





1 





2 


2 


1 


1 























-2 


























2 


2 





2 





2 


2 


2 





Total 


5 


2 


5 


2 


5 


4 


6 


2 



Table 6: The number of free coefficients in the j = 2 sector for liomogeneous 
turbulence and for different geometries 



Now we wish to use this form to fit for the scaling exponent Q in the 
structure function S^^{r) from data set I where a = /3 = 3 and the azimuthal 
angle of r in the geometry is — 7r/2.: 



S^i2,k=i{r, e,^ = 7T/2) - a9,2,i^^^" [-2^'' (Cf ^ + 2) sin 6 cos 6 

+ 2Cf^(Cf^ -2)cos3 6lsin6'] 
5'i2,fc=-2('-,^,0-V2) = 

^'i2,fc=2('-,^,0 = V2) = la9,2,2r^^"[-2Cf (Cf^ -2)cos2 0sin2 0] 

+ ai,2,2r^^"[-2Cf'(Cf^ -2)sin2 0]. (59) 

We see that choosing a particular geometry eliminates certain tensor contribu- 
tions. In the case of set I we are left with 3 independent coefficients for m ^ 0, 
the 2 coefficients from the m = contribution (Eq. |5^), and the single coeffi- 
cient from the isotropic sector giving a total of 6 fit parameters. The general 
forms in ^ can be used along with the fc = (axisymmetric) contribution ^ 
to fit to any second order tensor object. For convenience, Table ^ shows the 
number of independent coefficients that a few different experimental geometries 
we have will allow in the j = 2 sector. It must be kept in mind that these forms 
are to be used only when homogeneity is known to exist. If there is inhomo- 
geneity, the incompressibility condition cannot be used to provide constraints in 
the various parity and symmetry sectors, and we must in general mix different 
parity objects using only the geometry of the experiment itself to eliminate any 
terms. 



B The j=l component in the inhomogeneous 
case 
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B.l Antisymmetric contribution 

We consider the tensor 

T"^(r) =< u«(x + r) - u«(x))(u^(x + r) + u''(x)) > . (60) 

This object is trivially zero for a = 0. In our experimental setup, we measure 
at points separated in the shear direction and therefore have inhomogcncity 
which makes the object of mixed parity and symmetry. We cannot apply the 
incompressibility condition in same parity/symmetry sectors as before to pro- 
vide constraints. We must in general use all 7 irreducible tensor forms. This 
would mean fitting for 7 x 3 = 21 independent coefficients plus 1 exponent ^2^' 
in the anisotropic sector, together with 2 coefficients in the isotropic sector. In 
order to pare down the number of parameters we are fitting for, we look at the 
antisymmetric part of T"'^(r), 

T«^(r) = ^-L^ ^-1 = (u«(x)u'5(x + r)) - (u''(x)u«(x + r)), (61) 

which will only have contributions from the antisymmetric j = 1 basis tensors. 
These are: 

• Antisymmetric, odd parity 

^3,1™ = r-'[r"d^ - r'5a«]rYi,™(f ) (62) 

• Antisymmetric, even parity 

Bti,^ = v-h^^'^d.rY.^^ir). (63) 

As with the j = 2 case we form a real basis rYi,fe(f) from the (in general) 
complex rYi^rni^) in order to obtain real coefficients in our fits: 

rYi,fe=o(f) = rFi,o(f) = rcos6l = ra 

~ yi,i(f) + Fi,i(f) 

Zi 

= r sin sin 4> = 

~ yi,-i(f)-yi,i(i^) 

ryi,fe=_i(r) = r— ■ 

= r sin ^ cos (j) = r\ 

(64) 

And the final forms are 
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53,ll(f) 


= r~ 






= r~ 


-2e"'3^r^(r.p) 


52,li(f) 


= r~ 






= r~ 






~ r~ 


-2e"^''r^(r.m) 




= r~ 





Wc now have 9 independent terms and cannot apply incompressibility in 
order to reduce them further. Instead, we use the geometrical constraints of the 
experiment to do this: 

• (j) = Q (vertical separation), a = 3, /3 = 1 

B|\o(r,6',(;A = 0) = -sin6' 

Bl\Ar,e,c^ = Q) = l 

Bl\_^^{r,e,(j) = Q)=cose. (66) 

There are no contributions from the reflection-symmetric terms in the j = 
isotropic sector since these are symmetric in the indices. The helicity term in 
j = also does not contribute because of the geometry. So, to lowest order, we 
have 

= aa^i^o?"^ (- sin 6*) + 02,1, ir^2 + as, 1,-1?'^ cos0. (67) 
As always, we have written the scale dependent prefactors aqjm{r) as having a 

power law dependence aqjmr'^^ ■ We have 3 unknown independent coefficients 
and 1 unknown exponent to fit for in the data. 

B.2 Symmetric contribution 

We consider the structure function 

5'"/3(r) (u"(x + r) - w"(x))(u'3(x + r) - u^ix)) > (68) 

in the case where we have homogeneous flow. This object is symmetric in the 
indices by construction. It is easily seen that homogeneity implies even parity 
in r: 

S'^^i-r) = S"^^{r). (69) 

We reason that this object cannot exhibit a j = 1 contribution from the 50(3) 
representation in the following. Homogeneity allows us to use the incompress- 
ibility condition 

d^S"'^ = 

d^S"/^ = (70) 
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separately on the basis tensors of a given parity and symmetry in order to 
give relationships between their coefficients. For even parity, symmetric case 
we have for general j > 2 just two basis tensors and they must occur in some 
linear combination with incomprcssibility providing a constraint between the 
two coefficients. However, for j = 1 we only have one such tensor in the even 
parity, symmetric group. Therefore, by incomprcssibility, its coefficient must 
vanish. Consequently, we cannot have a 7 = 1 contribution for the even parity 
(homogeneous), symmetric structure function. 

Now, we consider the case of an experiment when r has some component in 
the inhomogcncous direction. Now, it is no longer true that S'"'^(r) is of even 
parity. Moreover it is not possible to use incomprcssibility as above to exclude 
the existence of a j = 1 contribution. We must look at all j = 1 basis tensors 
that are symmetric, but not confined to even parity. These are: 
• Odd parity, symmetric 



-^5 1 fc (^) ^^'^ (f ) vanish because of the taking of the double derivative 

of an object of single power in r. We thus have 4 different contributions to 
symmetric j = 1 and each of these is of 3 dimensions {k = —1,0, 1) giving in 
general 12 terms in all: 




r5"5^rfifc(f) =0 



(71) 



• Even parity, symmetric 





^5"''(r-n) 

^r"r^(r • n) 

^[(r"m'5+r'5m")(r 

M"^(r-p) 

V"r'^(r -p) 
2[(r"m^ + r'^m")(r 
M"'^(r-m) 



n) 



P) 



(rV + rV)(r -m)] 



(r"n^ + r'^n")(r -m)] 
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k 


= O,a = /3 = 3 


= 0,a = /3 = 1 


= O,a = 3,/3= 1 


6* 7^ 






= 




6' = 





3 


3 


2 


1 


2 





1 


1 





1 











-1 


2 





3 





2 


1 


Total 


6 


3 


6 


1 


4 


1 



Table 7: The number of free coefficients in the symmetric j = 1 sector for 
inhomogeneous turbulence and for different geometries. 



a/3 



a/3 



r'^^{r ■ m) 

= r-2[(rV + rV)(r ' n) - (r°n'' + r''n")(r • p)]. 



(73) 



These are all the possible j = 1 contributions to the symmetric, mixed parity 
(inhomogeneous) structure function. 

For our experimental setup II, we want to analyze the inhomogeneous struc- 
ture function in the case a = [3 = 3, and azimuthal angle (j) = (which cor- 
responds to vertical separation). For that case, we obtain the basis tensors to 
be 



33 



33 



r33 
r33 



= COS0 

= 2cos6' 

= cos^ e 

= — 2 cos 6* sin ( 

= sin 

= cos^ 9 sin 0. 



(74) 



Table | gives the number of free coefficients in the symmetric j = 1 sec- 
tor in the fit to the inhomogeneous structure function for various geometrical 
configurations. 



C Tests of the robustness of the interpolation 
formula 

In order to test the robustness of the interpolation formula Eq. we performed 
the following additional calculations. We considered the data from the probe 
at the height of 0.54 m. For each order n of the structure function, we defined 
a 'window' of data extending over two decades of the separation scale, r. We 
first placed the lower edge of the window well inside the dissipation range and 
fit the interpolation formula to the data in the first window. We then moved 
the lower edge of the window by half a decade and fit the formula to the data 
in the next window. In this manner, we proceeded until the upper edge of the 
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last window corresponded to the largest value of r. The entire range of r yields 
five windows. We thus obtained five values of the parameter C„ and calculate 

(2) 

the scaling exponent Ck = n — 2C„ in each case, giving some indication of the 
robustness of our result. 



C2 


0.35 ±0.1 


0.35 ±0.05 


0.39 ±0.02 


0.38 ±0.05 


0.38 ±0.07 


S2 


1.31 ±0.2 


1.30 ±0.10 


1.21 ±0.04 


1.24 ±0.10 


1.23 ±0.14 




Table 8: Second order: ^ 




= 1.25 ±0.05 






C3 


0.99 ±0.03 


0.95 ±0.04 


0.88 ±0.07 


0.91 ±0.04 


0.96 ±0.08 


A2) 
^3 


1.01 ±0.06 


1.10 ±0.08 


1.3 ±0.14 


1.2 ±0.08 


1.1±0.16 






Table 9: Third order: 




1.14±0.11. 






C4 


1.21 ±0.07 


1.12 ±0.09 


1.15 ±0.03 


1.29 ±0.1 


1.21 ±0.08 


S4 


1.58 ±0.14 


1.76 ±0.18 


1.7 ±0.06 


1.42 ±0.2 


1.58 ±0.16 



Table 10: Fourth order: ^ = 1.61 ± 0.13. 



Tables p[jl0| present the results of performing these checks on structure func- 
tions of the second, third and fourth order. The mean and standard deviation 
of the exponent values are given in the caption for each table. It is found that 
the mean value in each case is in close agreement to the value of the exponents 
presented in the main text which were obtained by a fit to the entire range of 
data. This gives us greater confidence in the use of the interpolation formula. 
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